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Abstract 

We systematically compute the effective short-range potentials arising from second order QCD- 
diagrams related to bound states of quarks, antiquarks, and gluons. Our formalism relies on the 
assumption that the exchanged gluons are massless, while the constituent gluons as well as the 



q . lightest quarks acquire a nonvanishing constituent mass because of confinement. We recover the 

well-known Fermi-Breit interaction in the heavy quark limit. Such an effective potential has been 
proved in the past to be relevant for the building of accurate potential models describing usual 
hadrons (mesons and baryons). The general potentials we obtain in this work are also expected to 
be useful in the understanding of exotic hadrons like glueballs and qqg hybrids. In particular, we 
compute for the first time an effective short-range quark-gluon potential, and show the existence 
of a quadrupolar interaction term in this case. We also discuss the influence of a possible nonzero 



mass for the exchanged gluons. 
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I. INTRODUCTION 



Mesons and baryons have been intensively studied in the framework of potential models 
for a long time l|, |2|. Moreover, since QCD is a nonabelian gauge theory, the gluons are able 
to form a bound state, at least in principle. These bound states of two, three or more gluons 
are called glueballs. The theoretical computation of the glueball spectrum is a problem 
which currently deserves much interest, either in lattice QCD j^], or within the framework 



of effective models 



Other so-called "exotic" hadrons are hybrid mesons: 



Mesons in which the color field is in an excited state. They can be seen as a three-body qqg 
state, with a constituent gluon. Again, hybrid mesons are the subject of many works, in 
lattice QCD 10, 11]. and effective models 12, 13, Q, Q, l^]. The experimental detection 



of exotic hadrons is also an active field of research, although no definitive conclusion has 
been drawn |17|. The confrontation between theoretical approaches and experimental data 
appears thus particularly interesting in order to confirm or invalidate some exotic candidates. 

Let us now turn our attention to effective QCD models. It is a well-known fact that 
the short-range interaction between two heavy quarks is given by the Fermi-Breit poten- 
tial IB, P- 395]. This corresponds to the one-gluon exchange process between the quarks, 
and contains fine structure terms like spin-orbit interaction and tensor force for example. 
Effective models with a linear confining potential and the Fermi-Breit potential have been 
showed to correctly describe the meson and baryon spectrum for a long time [1]. Apart 
from phenomenological arguments, the static potential between a quark and an antiquark 
has also been intensively studied in numerous theoretical works. At large distances, typically 
for distances larger than 2.5 GeV -1 , the static potential indeed appears to be linearly rising 
(see for example the lattice calculations of ref. [3]). At short distances however, the static 
potential has been shown to be no longer linear and nonperturbative but rather perturbative 
and Coulombic, just as is the static limit of the Fermi-Breit interaction |20j. 

Logically, the same kind of models - linear confinement and Coulomb term at long- and 
short-range respectively - has since been applied to describe exotic hadrons, like glueballs 
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9] and hybrids [12], ]13|, iJJJ. However, if the derivation of the Fermi-Breit 



q, y, y, 0, y. 



potential between two quarks is a standard procedure, it is not the case for the interactions 



between two gluons (see for example the discrepancies between refs. 



). Moreover, 



to our knowledge, the equivalent of the Fermi-Breit potential between a quark and a gluon 
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has not yet been computed. Since recent references, based on lattice QCD computations, 
validate the picture of a glueball as a bound state of constituent gluons [2l| and of a hybrid 



meson as a qqg system JjJ, the accurate knowledge of the gluon-gluon and quark-gluon 
potentials appears to be crucial point to further understand the properties of these exotic 
hadrons. It is thus of interest to compute all these potentials in a systematic way, within a 
single formalism, in order to obtain a set of potentials which should be applicable to most 
hadrons, either exotic or not. These computations are performed in this work. 

Our paper is organized as follows. In sec. [Til we introduce the general procedure for 
extracting the effective potential corresponding to a particular matrix element. Before ap- 
plying it to all the possible second order QCD-diagrams related to bound states, we have to 
know the non relativistic wave functions of the in and out states. This is presented in sec. IIHI 
Then, we compute the effective potentials for the diagrams involving pairs of quarks and 
antiquarks in sec. ITVl In sec. |V] we consider the interactions between two gluons, and the 



uons 

,3), 



interactions between (anti) quarks and gluons are treated in sec. I VI I The exchanged 
are massless in our formalism. Since it is not the case in some approaches (see refs. [4], 
we then discuss in sec. I VIII the influence of a possible nonzero mass for the exchanged gluons. 
In sec. IVHH we discuss the validity of the potentials which are computed. Since we want our 
potentials to be relevant in various systems - exotic or non exotic hadrons -, we compute 
in sec. HXI the color factors appearing in the various potentials we obtain. Finally, we draw 
some conclusions in sec. Kl We also added two appendices, where we recall the Feynman 
rules for QCD (appendix |A]) and some useful Fourier transforms (appendix IB1. 

II. GENERAL METHOD 

Let us consider a general Feynman diagram of order g 2 , with g the strong coupling 
constant 

1,A I' A' 




2,B 2'B' 



A,B,A',B' denote the color indices of the particles 1, 2, 1', and 2'. As the goal of this 
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work is to compute effective potentials for bound states of quarks, antiquarks, or gluons, we 
will restrict ourselves to diagrams where the initial and final particle content are the same: 
1 = 1' and 2 = 2'. The application of Feynman rules in momentum space, summed up in 
Appendix |Al leads us formally to the following matrix element (in units where h = 1) 

M fi = -z(2vr) 4 5 4 (P - P') M(*f a ; g), (1) 

where q is the exchanged 4-momentum, and ^p a , ^o'b are the properly normalized wave 
functions of the in and out states. In our notations, a, (3 = 1, 2 label the particles, C = A,B, 
and C = A', B' . These wave functions can be schematically decomposed as 

^ = c ®^(p), (2) 

with cj) c the color wave function, and ip(p) the spatial wave function in momentum space, 
p being the 4-momentum of the particle. The explicit form of these wave functions in the 
nonrelativistic limit will be given in sec. IHIl Then, the crucial and most complex step of 
the procedure is to rewrite M as 

M = «> [g 2 O A , B , AB U(q , q,p a , S a )] (3) 

p a and S a are the momentum and spin of the particle a. In the latter, bold quantities 
will always denote a vector. The spin operators appear through the rearrangement of the 
wave functions, and Oa'B'Ab is a color operator emerging from the vertices. We will show 
in sec. IIXI that this operator turns out to be a real number when correctly contracted with 
the color functions of the particles, <p c . We denote C this color factor. 

Once the decomposition is known, the effective potential we want to find is readily 
obtained in position space by a Fourier transform, denoted as J 7 , which changes q in r, the 
relative position of the two interacting particles: 

V(r) = Ca s T [4nU(q ,q, P p,Sp)} . (4) 

We used the definition g 2 = 4iras, which has the advantage of simplifying the Fourier 
transforms. All the Fourier transforms we need are listed in Appendix [B] 

We have at last to make a remark which will be valid in all our work: The effective 
potentials, and all the other quantities appearing in the latter, will be computed at the 
order c~ 2 , with c the speed of light. This is a standard approximation scheme, justified 
because the most relevant contributions as spin-orbit, tensor interaction, . . . , arise from 
these first relativistic corrections. 
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III. IN AND OUT STATES 



A. Quarks and antiquarks 

Quarks and antiquarks are spin- 1/2 particles, so they are described by Dirac spinors, 
denoted as u l (p) for the quarks and v^(p) for the antiquarks, with p the momentum of the 
particle and i,j their color indices. The wave function of a Dirac spinor, whose mass is m, 
can be computed at the order c~ 2 thanks to a Foldy-Wouthuysen transformation 22j. In 
our notations, it can be recast in the form (J2J), and reads 

u\p) = 0*®u(p), 

v j (p) = <f) j ® v (p), (5) 

with <b k the color functions, and 



u(p) = V2m 



[1 — p 2 /8m 2 c 2 ) w 
(cr • p/2mc) w 



(cr-p/2mc)w 

v(p) = V2m \ | , (6) 

'1 — p 2 /8m 2 c 2 ) w 

where w is a two component spinor and cr are the Pauli matrices (see Appendix [X]). In order 
to rewrite the matrix element on the form ([3]), we will have to isolate a factor \f2mw for 
each wave function. 



The following useful relations hold for both u(p) and v (p) |18|, p. 281] 

U( P V«(P) = 2mu>» (l - g2 + 2 ^ 2c ( 2 gXP) ) », (7a) 

i\ ± f ier x q q + 2p\ 

u(p ) 1U {p) = 2m u,t + j W , (7b) 

with q = p' — p, u = vfj . 

An important relation between the Pauli matrices is 

(crA)(crB) = A B + icr ■ (A x B), (8) 

in which A and B are two vectorial operators. The spin operator of a quark or an antiquark 
is defined by S = a/2. 
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B. Gluons 



A massive, colored, spin-1 particle is described by the wave function 

E a v (p)=<p a ®E u ( P ), (9) 

where e M is a polarisation vector satisfying the Proca equation with the transverse condition 

p v e v {p) = 0. (10) 

Relation ffTUj) together with the normalisation condition e u e u = — 1 leads to the nonrelativis- 
tic expression 

V mc 2m l c l / 

with 

p = ( mc + -?—,p) , and e 2 = 1, (12) 

in agreement with refs. |5j, |7| • The polarization vector e does not depend on the momentum 
p anymore. The factor to isolate in the matrix element will now be given by ^/2p° e, by 
analogy with the Dirac spinors. The spin operators for a spin-1 particle are defined by the 
matrices S, where 

£ kim j g ^ ne L ev i_Qi v j^ a symbol. The following useful relation involves the spin-1 matrices 

A t B, = [A ■ B - (S ■ B)(S ■ A)] tJ . (14) 

For QCD applications as those we consider here, the spin-1 particles which have to be 
taken into account are the gluons. These are massless particles, so, one can wonder if formula 
( TTTT) makes sense. In our framework however, exchanged gluons, which are true massless 
gauge particles, have to be distinguished from the constituent gluons of the in and out states. 
These ones are confined inside a hadron, and one can consider that even if the bare gluon 
mass is zero, they acquire a nonzero constituent mass due to the confining interaction, m 
has thus here to be understood as the constituent gluon mass in a particular bound state, 



ras a value around 0.5 — 0.6 GeV 
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231 ] . This argument has already been 



which typically 

used in refs. p, Il2| . and will be further developed in sec. IIV CI Let us note that the same 
argument holds in the case of the massless u and d quarks, whose constituent mass is around 
0.3 GeV Q. 



IV. INTERACTIONS BETWEEN QUARKS AND ANTIQUARKS 
A. Fermi-Breit potential 

As an illustration of the general method presented in the previous section, we begin by 
the best-known diagram, that is the one-gluon exchange between two quarks. Applying 
the Feynman rules of Appendix we can write the matrix element corresponding to the 
diagram 



P 




P2,k v,b p' 2 ,l 



in the form (pQ), which is 

M fi = -i(2vf5\P' - P)(-/)C» 1; ^ J{\q)D^(-q). (15) 
The color operator is given by 

0W« = f ^ (16) 
and the conserved current of the quarks by 

Jtt(q)=u' a (p + q)^u a (p). (17) 

We omitted to note the color functions in order to simplify the notations. The explicit 
computation of the color factors will be done in sec. IIXI Let us call C\ the color factor 
corresponding to Oi^i- 

As the current J^io) is obviously known thanks to the relations (JTj), we only need to 
specify the form of the propagator D^ v in order to develop eq. (fT5l) . Let us note that 

Q=P2-P2=Pl-Pl ( 18 ) 

An interesting gauge choice is then to set 9 = in the propagator (lAlj) . This leads to the 
expressions 

A)o = -^, A)* = 0, D ik =±L k -^f\. (19) 
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Putting all these elements together, one can rewrite eq. (fT5"j) in the form ©, i. e. 

M= (im 1 m 2 )w'^ [g 2 C 1 U(q, Pl ,p 2 )] Wl w 2 , (20) 



with 



U{q,Pi,p 2 ) 



2 n 2 



Sm\c 



2 n 2 



?cri • (q x pi) «<x 2 -{qxp 2 



4m\c 2 q 2 



+ 



(gl ■ g)(P2 • g) 

mim 2 q 4 c 2 
(cti • qf)(<r 2 • q) <r, ■ a. 



Pi • p 2 uti ■ [q x p 2 



Arrive 2 q 2 

icr 2 ■ (q x pi) 



m\m 2 q 2 c 2 



2m,\m 2 c 2 q 2 



2m\m 2 (?q 2 



(21) 



Amim 2 c 2 q 2 Amim 2 c 2 
In the latter, the c factors will not be written anymore in order to simplify the notations. 
The last step is to compute the Fourier transform of U, as formulated in eq. (j3J). By using 
the relations of Appendix [Bj we obtain the following effective potential 



V{r) = da s 



n 












\m\ 


mi J 



2m 2 r 3 2m 2 r 3 



1 



2m\m 2 r 
1 



P\ P2 



(Pi • r )(P2 • r) 



1 



+ 



m 1 m 2 r° 



Si ■ S 2 — 3 



[Sr-rXSz-r) 



m 1 m 2 r J 
8tt 



3mim 2 



\L\ ■ S 2 — L 2 ■ Si] 



Si ■ S 2 8 3 (r) 



(22) 
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, p. 395]. It is 



where Li = r x p^. This is the well-known Fermi-Breit interaction 
interesting to rewrite this contribution by using 

l = ~ mA , n = £ 1 + £ 2 , 

mi + m 2 

which are simply the relative and center of mass variables of the interacting pair, and also 
by introducing new spin variables as follows 



(23a) 



S = S\ + S 2 , 6 = Si — S 2 



(23b) 



Potential (1221 then becomes 



V(r) = Cia s - 



71 



6 -AS 2 



1 



\ 3mim 2 



m\ + m\ + 4mim 2 



4m\m\ 
[m 2 - mi) 



2m 2 
L S 



2m\ 



+ 



5 3 {r) 



2mim 2 



P1P2 , (pi • r)(pa ■ r) 



Amim 2 r 3 



S-Q 



4(mi + m 2 )mim 2 



+ 



2mim 2 



+ 



(S 



Am^mf J r 



L 6 



(24) 
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It can be checked that the skewsymmetric spin-orbit interaction vanishes for two particles 
of equal masses. Indeed, in this case, the wave function must either have a good symmetry 
under the exchange 1 •*->• 2 if both particles are identical, or a good C-parity for a particle- 
antiparticle pair. In a two-body system, O = can be imposed by working in the center of 
mass frame. 

If we considered an antiquark interacting with a quark or an antiquark, we would have 
found the same expression (|22p . but with a different color factor, since the color operator 
would be different. These possible color operators are listed in eq. ( 187)) . and denoted as 
2 ,3-,ijki respectively. 



B. Annihilation diagram 

In the particular case of the interaction between a quark and an antiquark of the same 
flavor, the contribution of a new diagram must be added, that is the annihilation diagram: 




The flavor of the outcoming pair is assumed to be identical to the incoming one. 
With the same indices as for the scattering diagram, we find 

Mfi = - 2(2vr)V(P / - P){-g 2 )0^ kl 

(25) 

x [v 2 ^ mi] £V [u[Y v' 2 ] , 

with 

O^-ijki = ^(A a )jfcj(A a )j7. (26) 

The independence of the matrix element ( 125)) in 9 is readily checked. Consequently, the 
relevant part of the gluon propagator in the s-channel is given by 

v ~4m 2 ' 1 j 

because we can suppose that the momentum of the created gluon is very smaller that its 
rest energy. The wave functions (E]) can be taken at the lowest order, that is 




u = V2m I I , v = V2m I ° I . (28) 



w 



After a rearrangement of formula (1251) thanks to the relation 



1 3 X 



one finds 



M = (w'£) C 4 -^ [3 + o-i ■ cr 2 ] (w 2 ) c wi- 



(29) 



(30) 



O^ijki has been replaced by the corresponding factor C 4 , and we defined the charge conju- 
gation as 

w c a = ia 2 w* a . (31) 
After a Fourier transform, we get the annihilation potential 



v(r) = a 



27TCts 



4 



+ Si ■ S 2 



<S 3 (r), 



(32) 



which has logically the form of a contact interaction. Moreover, it is a projector on a spin-1 
state: The quark-antiquark pair must have the quantum numbers of a gluon in order to 
make the annihilation possible. The color factor C 4 is indeed also a projector, on the octet 
state (see table [III)- Let us note that our potential ( 1321) is in agreement with the one of 
ref. [24]. 



C. Comparison with background perturbation theory 



The spin-dependent corrections coming from QCD are computed with the background 
perturbation theory in ref. [3] for the qq pair in a meson. In our notations, the result reads 

Si ■ L\ S2 • L2 



V SD (r) 



2/x? 
Si ■ S 2 



VJr) 



2f4 J L 

1 



-<9 r e(r) + -<9 r T4(r) 



S 2 ■ Li — Si ■ L 2 1 , . 

H 9 r V2(r) 

T 



3 (gx-r)(g a -r)_ 



(33) 



Vs(r) 



3yUl/i 2 3yUlyU 2 

where e(r) and V^(r) are rather complicated functions of the quark correlators. Fortunately, 
at the order g 2 they take simpler expressions, which are 



V l (r) = 0, V 2 (r) = e(r) = -§sa 



4a,c 



V4H 



y7ra s 5 3 (r) 



(34) 
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One can check that, by injecting the relations (I34p into the potential (I33p . one formally 
obtains the spin-dependent part of the Fermi-Breit potential (1221) for a meson, that is for a 
color factor equal to C2 = —4/3. The only difference is that, in the background field theory, 
the bare masses m a are systematically replaced by fi a , which are interpreted as dynamical, 
or constituent quark masses. 

If one sees the potentials we are computing as perturbations of an effective Hamiltonian 
(typically a spinless Salpeter one with linear confining potential), the \i a can be defined 



as 



26] 



^ = (VP« + <), (35) 

where the averages are computed with the unperturbed wave functions. The advantage of 
this approach is that the [i a are always nonzero; the development of the matrix elements 
in powers of l//i Q is thus relevant, even for light particles. Consequently, the equivalence 
between potentials (122j) and (l3"3"l) . suggest that in all our effective potentials, the quark 
masses should be seen as the dynamical ones, the /i Q defined by eq. (l33|) . 

The same arguments also hold for gluons: Even if the exchanged gluons are seen as 
massless gauge particles, it is usually considered that the constituent gluons, in a glueball 
for example, acquire a mass because of the confining interaction. This point of view can 



also be checked with the background perturbation theory (see ref. js]), and has already 
given good results in the description of glueballs 8, sj. Thus, in our case, the "gluon mass" 
appearing in eq. ( JTTl) has to be understood as the constituent one. 



V. GLUON-GLUON INTERACTIONS 
A. Scattering 

The computation of the effective potential coming from the scattering diagram between 
two gluons, 
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is rather similar to the one of the scattering diagram between two quarks. The analog of 
formula (jT5j) is now 



M fi = -^2vr) 4 5 4 (P' - P)(-g 2 )0 5 , abcd G*{q)D^G v 2 {-q), 



(36) 



with the color operator 



and the conserved current 



5:abcd 



f baft 



edci 



(37) 



(38) 



With the aid of the definitions (ITT]) and (jT3"j) . and the property 



(a ■ e a )(b ■ e*) - (b ■ e a )(a ■ e*) = i(b x a) ■ S a (e a ■ e* a 



G a (q) = 2p° a 



-(2p a -q)-e a e a + e n e„ 



(39) 



we can write the different components of as 

G° a (q) = -2p° a e* |l - [i(q x p a ) -S a + q 2 - (S a ■ qf] } e a (40) 



2mc mc mc 

As the gluon propagator is given by relations (1191) . the matrix element f)36p becomes 



(41) 



M 



9' 



fi 



G1G2 — G\ ■ G2 + 



(q-G 1 )(q-G 2 ) 
q 2 



(42) 



where we used the obvious notation G\ = Gi(q) and G2 = G 2 (—q). By simplification of 
eq. (j4"2l . and with formula (THj) . one can arrive at the analog of eq. ( !20|) . that is 



M=(Ap lP 2 )ele; [g 2 C 5 U(q, Pl ,p 2 )] e 1 e 2 , 



(43) 



with 

U(q,Pi,P2) 



I I iSx ■ (q x Pl ) iS 2 ■ (q x p 2 ) (pj ■ qr)(p 2 • 9) Pi ■ P2 



2m 2 q 2 



2m 2 q 2 



m 2 q 4 



iS 1 -(qxp 2 ) iS 2 -{qxpi) {Si ■ q){S 2 • q) Si ■ S 2 

00 o o 



m 2 q 2 



m 2 q 2 



m 2 q 2 



(S 1 ■ q) 2 + (5 2 • qf 
2m 2 q 2 



m 2 q 2 



(44) 
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The constituent gluons are assumed to have the same effective mass: This is clear because 
of the symmetry of the problem. The effective potential in position space is thus 



V(r) 



C 5 <y s { - 

r 



47r r3/ \ Li ■ Si L 2 ■ S 2 
o (r) 1 — 

3m 2 2m 2 r 3 2m 2 r 3 

S\ ■ S 2 



1 



Pi ■ P2 + 



(Pi ■ r )(P2 ■ r) 



m 2 r 3 



1 [Li .S 2 -L 2 -Si] + -^ 



- 3 



2m 2 r 
{Si-r){S 2 -r) 



S7T 



3m 2 



Si ■ S 2 5 3 (r) 



m 2 r 3 



[(Si ■ rf + (S 2 ■ rf] } 



(45) 



2m 2 r 5 

This potential is very similar to ( 1221) . except in two terms. The first one is the contact 
interaction, which is stronger than in the quark-quark case, and the second one is, at first 
sight, a new interaction term, written in the last line of eq. (1451) . The existence of this term 
has already been pointed out in ref. j^], but here we go further than this reference because 
we completely remove the e a from our expression. Actually, these terms are a part of the 
tensor force, because if we rewrite the potential f j4"5l) with the new variables defined by ( f23l ). 
we obtain 



V(r) = C 5 a s 



' +i^(3-5 2 ) 5 3 (r) 



1 



r 3m 2 
3 L ■ S 1 6 ■ n 



2m 2 



Pi ■ P2 (Pi ■ r){p 2 ■ r) 



1 



2m 2 r 3 4m 2 r 3 



2m 2 



-=--3 



[S-rf 



(46) 



B. Annihilation 



The annihilation of two gluons, represented by the diagram 




is a special case because its contribution is always zero. First of all, the color operator 
associated to this diagram is 

G -abcd = fcafebd- (47) 

It is shown in sec. IIX B I that the corresponding color factor is non zero only if the gluon pair 
is in an antisymmetric color octet. But even in this case, a calculation similar to the one 
performed in sec. IIV Bl shows that the matrix element is zero at the order c~ 2 . Indeed, as 
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the gluon propagator in the s-channel can be taken of the form (1271) . the wave functions and 
momenta will simply be given by their lowest order approximation 

e = (0,e), p=(m,0), (48) 

respectively This approximation causes the matrix element to vanish in every case, as it is 
mentioned in ref. [7]. 



C. Seagull 

The seagull diagram possesses a unique vertex: The four-gluon one. 




Applying the corresponding Feynman rules together with the Jacobi identity 

fabfecd + facfedb + fadfebc = 0, (49) 

we obtain the following matrix element 

M fi = -*(2vr) 4 5 4 (P' - P)^- {0 6 , abcd {2t\ ■e 1 e* 2 -e 2 -e* 1 - e 2 e* • e 1 - e\ ■ e* e x ■ e 2 ) 

+ 5 . abcd (el ■ ei e* 2 ■ e 2 - e* 2 ■ e x • e 2 )} . (50) 

An additional 1/2 factor is present in eq. ( |50l) . This is due to the degeneracy of the different 
diagrams. Indeed, if we attribute a weight 1 to the gluon-gluon scattering diagram, we have 
to attribute a weight 1/2 to the annihilation and to the seagull diagrams j^]. 

The effective potential corresponding to the seagull diagram should clearly be a contact- 
like interaction. It is thus sufficient to inject the lowest-order expressions (l4"8"j) to simplify 
eq. ( |50l) . The following relations have also to be used 

(e* • e 2 )(e* 2 ■ ei) = e* 2 e\ (Si ■ S 2 ) e 2 e x + (e* • e 2 )(ei • e 2 ), (51a) 
Tr[(e;S Q ).(S Q e Q )] = e* a S 2 a e a . (51b) 

Moreover, 

(elS^SUelSfie.Sy = e^S, ■ S 2 ) 2 e 1 e 2 w 0, (51c) 
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since such terms are higher order relativistic corrections. The effective potential in momen- 
tum space reads 



U(q) 



{Ce[S 1 -S 2 -4)+C 5 [S 1 -S 2 -2}}. 



Introducing the total spin of the system, we get in position space 



V(r) 



2m 2 



|c 6 




+ c 5 


>-< 


} 











S 3 (r). 



(52) 



(53) 



The complete potential to use to describe the interactions between two gluons is thus the 
sum of potentials ( H6|) and ( 1531 . 



VI. (ANTI) QUARK - GLUON INTERACTIONS 

The understanding of hybrid mesons currently deserves much interest in experimental as 
well as in theoretical physics. From a quark model point of view, hybrid mesons can be 
seen as a quark-antiquark-constituent gluon bound state. In the framework of a potential 
model, the dominant par t of the interaction is a confining term which is typically given by 
a\x q — x g \+a\xq — x g \ [271 ] . Then, short-range interactions should be added in order to build a 
spin-dependent model. The quark-antiquark potential has been computed in Sec. HV] in this 
case, the qq pair is in a color octet and couples to the constituent gluon in order for the total 
system to be in a color singlet. But, the short-range potentials between the constituent 
gluon and the quark or the antiquark have, up to our knowledge, never been computed. 
Consequently, the knowledge of (anti)quark-gluon interactions is a missing element in the 
building of accurate potential models of hybrid mesons. That is why we compute these 
interactions in the present section, applying the same procedure as the one used up to now. 



A. Scattering diagram 



There are three diagrams to consider when one deals with quark-gluon interactions 28, 
p. 196]. We begin with the following scattering diagram, 




> > % 

id v,d p' g ,k 
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In the following, the indices 1, 2 will be replaced by g, q in order to clearly distinguish between 
the symbols related to the gluon and to the quark respectively. 

Applying the Feynman rules, we can write the corresponding matrix element in form (CQ), 



M fi = -z(27r) 4 5 4 (P' - P)g 2 7lbcjk WMD^rJ-q), 



(54) 



where the color factor is 



O 



7;bcjk 



A* 

Y 



acb ) 



(55) 



and where J q (—q) and G^(q) were already computed in eqs. (!T7j) and (IHHj) . Using the 
propagator (ITTJ1) . the matrix element can rewritten as 



M 



fi 



G n J n + G g ■ J q — 



{q-G g ){q-J q 



9 1 



(56) 



After the same rearrangements as those which were done in sees. |IV]and|V] eq. (|56p becomes 



M = (Ap° g m q ) e*w« [g 2 C 7 U(q, p g , p q )] ew, 



(57) 



with 



U(q,p g ,p q ) = -r - 



1 1 \ iS g - (qx p g ) iS q ■ (qxp q 



2m 2 g 4m^ 



2m 2 q 2 



+ 



+ 



Pg ■ q)(Pg -q) Pg Pg , ^ 9 ■ (9 x Pg) ^ • X p ) 



m g m q q q 



mgin q q A 



+ 



m g m q q z 



m q m g q z 



, (S g -q)(S q -q) _ S q -S g (S g ■ q) 2 (S q ■ q 



2m 2 q 2 



2m 2 q 2 



(58) 



m g m q q' m g m q 

One can check that the quark part of potential ( l58l) is equal to the one of eq. (121D . and that 
the gluon part is the same as in eq. f)44p . which is rather coherent. The last term of this 
expression is actually not explicitly written in the quark-quark case. It is indeed eliminated 
thanks to the properties of the Pauli matrices through the relation 

(S q ■ q) 2 /q 2 = 1/4; 

consequently, it is only a constant term, which will be part of the contact interaction after a 
Fourier transform. In the gluon-gluon scattering, the gluonic term in (S g • q) 2 was naturally 
included in the tensor force. In the quark-gluon case, the masses m g ,m q are not equal, and 
this term has thus to be treated separately. 
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In position space, the effective potential is given by 
1 



V(r) = C 7 a s 



27T 71 

+ 



3m 2 2m 2 



6 6 (r) 



Lg ■ Sg , L q - S q 



■in 



'Im-f''' 2m 2 r 3 3m n m n 



2m g m q 



m g m q 



PgPg (Pg-r)jpg-r) 



1 



m g m q r 



Ug,,u q 



3 i L g ■ Sq~ L q - S g ] 



Sg-S q z (S g -r)(S q -r) 



3 (S g r) 2 



m 2 r 3 2m 2 g r 5 



(59) 



Rewriting it with the appropriate variables we have 



V(r) = C 7 a s 



+7T 



1 



r 2m g m q 
11 2 



Pg Pq , (Pg-r)(p g -r) 



3m g m q 3m 2 2m 2 3m g m q 



S 2 5 3 (r) 



m q + m g + Am g m q \ L ■ S 



4m 2 m 2 



1 d ■ CI (m q - nig) S ■ Q ( m 2 - L ■ S 



Am g m q r 3 Am q m g {m q + m g ) r 3 



2m g m q 



u q" u g\""q 

:s-r) 2 - + 



Am 2 m 2 / r 3 



m q - m g 
m 2 m q 



Q(r) 



with 



Q(r) 



1 3(S g ■ r) 



2r 5 



(60) 



(61) 



The physical meaning of this term can be clarified by using the following relation involving 
the spin-1 matrices {29), p. 55] 



(62) 



Qik are the polarisation matrices, defining the quadrupole tensor. Consequently, eq. ( 16TT) 
can be rewritten as 

Q(r) = (63) 



2 r 5 

A quadrupolar interaction is thus present in the potential. Such a term is logically absent 
in symmetric configurations, i.e. quark-quark and gluon-gluon systems, but appears in this 
asymmetric quark-gluon system. 

B. "Compton" diagrams 



The two remaining diagrams are analog to those of the Compton scattering in QED. We 
begin by the following annihilation diagram, 
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whose matrix element reads 



qz _ m z 



where the color factor is 



The quark propagator reads 



0%;bcjk — tC^*^ )kj- 



+ m) {m g + m q ) r y + m 1 



q 2 — m 2 



(64) 



(65) 



(66) 



(m g + m q ) 2 — m 2 ' 

in the approximation q = (m g + m q , 0). fh is the mass of the propagating quark, which we 
choose to be equal to m q . Using the lowest order approximations for the wave functions and 
the momenta, and thanks to formula 



(e' • cr)(e • a) = e' ■ e — e'(S ■ <r)e, 



(67) 



we find 



M = 4p° g m q w' ] e 
The effective potential, 

2nC 8 a s 



C s g 2 ' 2S,rS,! 



2m g (m g + 2m, 



u J 



we. 



V{r) 



15 
T 



-s 2 



S 3 (r) 



(68) 



(69) 



m g (m g + 2m q ) 

is a projector on a spin 1/2 state because the quark-gluon system must have the quantum 
numbers of a quark in order to form a single quark. If one now considers an antiquark 
instead of a quark, the only modification in the potential (169]) will be the color operator, 
given by 

Og-bcjk = -^{^ b ^ C )jk- (70) 

The second Compton-like diagram is 
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This diagram is no longer an annihilation one. It can indeed be seen as an exchange of a(n) 
(anti)quark between a (anti)quark-gluon pair. We will thus work in Breit's frame, as in the 
case of the usual scattering diagrams: q = (0,q). The exchanged quark has a constituent 
mass m q , which is always nonzero in our formalism. Consequently, we expect to find an 
expression in which a Yukawa potential appears, whose contribution will in every case be 
less important than the one of the quark-gluon scattering diagram since the exchanged 
gluons are massless. This leads us to assume that only the lowest order expression of the 
effective potential will be relevant. One obtains easily for this potential 

r 



V(r) = C 10 a s 



Cio corresponding to the color operator 



2m n 



S z - 



(71) 



(72) 



Let us remark that if one deals with a quark instead of an antiquark in this diagram, the 
only change in potential ( J7TT) will be again the color operator, given by 

1 



O 



XX,bcjk 



7(A c A 6 ) ifc . 



(73) 



VII. GLUONS: MASSIVE OR NOT MASSIVE 



In the computation of all the potentials that we presented up to now, the constituent 
gluons were considered as massive particles because of confinement, but the exchanged 
gluons remained massless gauge particles. The advantage of this assumption, already used 
in ref. jsj, is to treat the diagrams involving constituent (anti)quarks and/or gluons on an 
equal footing. Moreover, it allows to compute the potentials with the standard Feynman 
rules corresponding to the QCD Lagrangian. 

However, in some other approaches related to potential models of glueballs, it is argued 
that the exchanged gluons have to be massive as well in order for all the gluons to have an 
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equal status 0, Q]. In these references, consequently, a gluonic mass term has to be added 
to the QCD Lagrangian. For completeness, we discuss here the consequences of such a 
term. Actually, the Feynman rules will not be modified at the level of the gluon propagator: 
the Coulomb form in position space, denoted hereafter as U(r) = 1/r, will turn into the 
propagator of a massive vector particle, typically U(r) = e~ mr 'jr. Thus, at the first order in 
the relativistic corrections, only the scattering diagrams will be affected by this modification. 
The corresponding scattering potentials, denoted here as V s , and given by eqs. 
and (EDI), can be rewritten as 



V* = c ias { u 



/3-2S 2 



V 6mi?7i2 8m? 8m 2 



+ 



m\ + ra\ + 4mim2 
4m 2 m 2 > 



1 «"«2 

l-s u -- 

r 4mim2 



V 2 U 
1 



2mirri2 
r 



U' 

P\Up 2 - (pi ■ r)—(r-p 2 



(m 2 - mi) _ jy_ 
4(mi + m 2 )mim 2 t 



1 



6mim2 



U" 



r 



mi — m-2 
4m 2 mf 
JS-r) 



r 

L-S U - 

r 



(74) 



V, 



3 



3m 2 



(3 - S 2 ) V 2 U 



Ira 2 r 



2m 2 
1 



4m 2 r Qm 2 



P\Up 2 - 

u"- u - 

r 



( \ U '( \ 
(pi -r) — (r -P2) 



:s-r) 



, (75) 



Vi = C 7 a s {U- 
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2m g m q 



U' 

PiUp 2 - (pi ■ r) — (r-p 2 



11 



Ylm g m q 



6m 2 



Am g m q 



6m g m q 



-6 ■ n 



U' 



u" 



r 



3m g m q 



m 2 + m 2 + 4m„m, 



^2jj I " u q 1 ""g 1 "";/'"'/ 



4m 2 m 2 



4m 2 m 2 



m"„ \ U' 
L ■ d 

r 4(m 9 + m q )m g m q 

S-r) 2 



s n 



L 

U' 



2m 2 m q 



U" 



U'\ rQr 

r 



U' 



(76) 



The symmetrizations on the noncommuting operators have not been explicitly written in 
order to clarify to formula. It appears that all the scattering potentials we compute almost 
share the same structure, the only difference coming in the coefficients of V 2 t/, typically 
responsible for the contact interaction. The quadrupolar term is trivially zero except for the 
(anti)quark-gluon interactions. In this last case, we denoted the gluon and the (anti)quark 
masses by m g and m q respectively, in order to avoid confusion. It is worth mentioning that, 
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if Ay is the coefficient of V 2 £7 in the scattering potential Vu, it can be checked that 

A ^ ^ ' (77) 

qq 8m 2 8m 2 , 3m 1 m 2 ' 

A = * L _ 25 ^ 5 ^ (78) 

99 6m 2 6m 2 3m 2 ' 

A = - — - 2Sq ' Sg . (79) 

99 6m 2 8m 2 3m g m q 

It appears obviously that the differences only depend on the nature of the interacting parti- 
cles, each quark bringing a contribution — 1/Sm 2 , while this contribution is equal to — l/6m 2 
for a gluon. 

Clearly, the Coulomb form of U comes from the propagator of a massless particle. Fol- 
lowing ref. [3lJ, the expressions ( BP , ( 1751) . and (1751) are valid for any U(r), this function 
being the Fourier transform of the propagator of an exchanged vector particle. Thus, for 
a massive gluon, one should simply replace U(r) by a Yukawa potential U(r). Moreover, 
following eq. (lB7j) . V 2 U = m 2 U - 47r5 3 (r). 

Various gluon-gluon potentials involving massive gluons can be found in the literature (see 



Suon i 
,3), 



for example refs. [4, LDj), but they are neither equal between each other, nor equivalent to our 
potential. In recent works jsoj], the coefficients of the different terms (L • S, L • 8, ...) of the 
potential between two heavy quarks was found numerically by lattice QCD computations. 
These fitted coefficients agree with the Fermi-Breit potential. It could be interesting to 
find, by the same procedure, the same coefficients for two static sources in a color octet. 
This could lead to identify the most relevant gluon-gluon potential. Moreover, an accurate 
computation of various hadrons spectra will be able to check the relevance of one approach 
or another, through a comparison with experimental data or with lattice QCD. 



VIII. VALIDITY OF THE APPROACH 
A. The strong coupling constant 

The effective potentials we computed up to now involve Feynman diagrams of order g 2 , 
without loop. It is well known however that the inclusion of the one-loop diagrams leads 
to the conclusion that the strong coupling constant depends on the exchanged momentum 
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through 

(33 - 2N f ) In (£) 

where Nf is the number of quark flavors whose masses are lower than q 2 , and where A 
is the famous lambda-QCD parameter, whose value is around 300 MeV. Clearly, as{q 2 ) is 
smaller and smaller for increasing q 2 , QCD becoming perturbative a t hig h momentum. As an 
illustration, it is nowadays accepted that as(m 2 z ) = 0.1176±0.0020 [32|, with rriz ~ 92 GeV. 
For q 2 pa A 2 however, QCD is nonperturbative at low momentum since the strong coupling 
constant (IHUl) blows up. This low-momentum domain precisely concerns the bound states 
physics we are interested in. Consequently, it appears necessary to wonder whether the 
potentials we obtain are valid or not to describe hadronic bound states. 

Lattice QCD clearly shows that the static potential between a quark and an antiquark is 



nicely fitted by a Cornell shape, that is ar — 4a,s/3r, for a ph 0.20 GeV 2 and as ~ 0.22 33 



P- 



42]. Such a form suggests that the total energy between the quarks can be roughly separated 
into a confining (nonperturbative) part at long-range (ar), and a "residual", short-range, 
Coulombic part (— 4ag/3r). In this picture, as should be interpreted as a small, effec- 
tive, strong coupling constant. We point out that, although the Cornell potential correctly 
matches the data, it should only be seen as an interesting approximation of the real static 
potential, for example because it neglects the interplay between perturbative and nonper- 
turbative effects at intermediate distances. Moreover, the calculations of ref. [3] forbid 
the presence of a linear term at short distances. But, the Cornell potential is clearly the 
simplest way to interpolate between a Coulomb potential at short distances and a linear 
confining one at long distances, both asymptotic behaviors being confirmed by theoretical 



studies 
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20]. Anyway, the short-range part of the static potential appears to be mostly 



given by a perturbative Coulomb term, i.e. the static approximation of all the potential we 
computed up to now. Let us note that strong arguments indicate that this picture is also 
valid not only for bound states of heavy quarks and antiquarks, but also for bound states of 
massless gluons |2l|, thus in the light particle sector. 



Several attempts have been made in order to compute as(g 2 ) at all orders, especial 



solving the Dyson- Schwinger equations (see a rather complete list of references in ref. 34j ) 
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All these approaches qualitatively lead to 

as(q 2 ) = ^7 r > (81) 

(33-2AT / )ln(^))' 

where £,(q 2 ) is a monotonic function such that £ 2 (0) > and £ 2 (g 2 — > oo) — > 0. However, 
the explicit formula giving £(g 2 ) is different following the different works. Equation (1HT|) 
states that the strong coupling constant remains finite for q 2 = A 2 , and tends to a maximal 
value 015(0) - one speaks of "freezing" of the strong coupling constant. The value of 05(0) is 
still a matter of research, but no theoretical argument forbids 05(0) < 1, and some studies 
even allow for 015(0) to be as low as 0.4 34J], just around the typical values which are used 
in potential models 35]. Such a freezing scenario is coherent with the perturbative point 
of view that we adopted, and justifies in particular the assumption that the lowest order 
Feynman diagrams already give a relevant picture of the short-range interactions between 
two confined particles. 

For completeness, it should be pointed out that, although providing a coherent framework 
to understand to success of potential models and how to build them, the freezing of the strong 
coupling constant is not universally accepted, see for example the lattice computations of 
ref. 36] in which no freezing is observed. 



B. Nonrelativistic expansion 

Another important ingredient of our computations is a nonrelativistic expansion of the 
"wave functions" (Dirac spinors and polarization vectors) in powers of p 2 /fi 2 . Such an 
expansion is clearly valid for heavy quarks since (p 2 /fi 2 ) ~ (p 2 /m 2 ) <C 1. In this last case, 
our quark-quark potential reduces to the usual Fermi-Breit one, which is widely used in 
potential models. 

For light quarks and massless constituent gluons, it is readily checked that (p 2 //i 2 ) ~ 
1, and the nonrelativistic expansion is formally no longer valid. However, although no 
theoretical argument is able to justify it, the extension of the Fermi-Breit potential to light 



mesons and baryons has been £ 
with a quite good agreement 



lown to reproduce the experimental data for a long time 



0,3 



35]. Even the glueball spectrum, computed within a 
potential model where we used the gluon-gluon potential that we obtained in this work, is 
in agreement with lattice QCD when massless constituent gluons are used [351 ] . By doing 
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explicit computations, one is thus led to the empirical conclusion that the first relativistic 
corrections in p 2 /fi 2 contain most of the physically relevant information. But, it is clear 
that light quarks are at the limit of the validity range of potential approaches. 



For w 



of ref. {4], several studies have given theoretical arguments favoring massive constituent 
gluons, with a mass which could be as high as m g = 0.8 GeV. In this case, our gluon-gluon 
potentials remain valid, but with U(r) = e _msr /r, as we argued in Sec. IVHI At such high 
masses, (p 2 /ttv^) becomes small enough for the nonrelativistic expansion to be justified. 
The introduction of massive constituent gluons, in gluon-gluon as well as in (anti)quark- 
gluon potentials, could then bypass the conceptual problems caused by the nonrelativistic 
expansion. 

IX. COLOR FACTORS 

A. Clebsh-Gordan coefficients 

We have now to compute the value of the various color operators O^a'b'ab that we 
introduced in the previous sections. Up to now, we did not impose the two interacting 
particles to be in a particular representation of SU(3). Let us illustrate this with the qq 
interaction in a baryon. The color functions of the quarks were of course chosen to be in 
the fundamental representation of SU(3), but we did not explicitly asked that the qq pair 
had to be in the conjugate representation in order to make a singlet when coupled to the 
third quark. To include this point in our discussion, we have to contract the color operators 
with the appropriate Clebsh-Gordan coefficients (CGC) of the in and out states. Let A, B 
be the generic color indices of the two incoming particles, and A', B' those of the outcoming 
particles. In all our potentials, O^a'B'Ab has thus to be replaced by 

Cn = Fq A B O n -,A' B 1 AbFq B 1 (82) 

ABA' B' 

where F~ B are the CGC. C stands for the color index of the coupling, which is the same 
for the in and out states (the color is conserved during the interaction). In every cases, one 
can check that the result is independent of C . 

Let us now compute the F~ B , in a similar way as it is done it ref. 33, p. 51]. Let 



rat concerns the constituent gluons, it is worth mentioning that, since the work 



be the color wave function of a particle in the representation X of SU(3). In other 
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words, 5?(j) x = ie ■ T x (j) x with T x the generators of the representation X and e defining 
an infinitesimal gauge transformation (in this section, the vectors will be denoted with an 
arrow for convenience). Because the Fq are general CGC that couple the representations 
X and Y to form the representation Z, we can find them by requiring 

Si M^VZ] = ie- {f z ) c D [<f> x F^] . (83) 
Moreover, the CGC computed thanks to eq. fl83l) have to be normalised with the relation 

^ F *AB F AB = 6cD (g4) 
AB 

Our method can be in principle applied to any representation of SU(3). But in this work, 
we will restrict ourselves to systems involving only the representations 1, 3, 3 and 8, which 
are particularly relevant for hadron spectroscopy. In particular, the quarks, antiquarks and 
gluons are in the 3, 3, and 8 representations. Some CGC are computed in table [B As it 
is shown in the second column, they can be used to describe mesons, baryons, glueballs, 
hybrids, etc. We recall that e\ 3 is the Levi-Civita symbol, and that the constants f£ b and 
d ab are defined by the well-known relations 

[\\\ h \ =2iJ2f?* C , (85) 

c 

{A fl ,A 6 }=2VrffA c + V. (86) 



B. Computation of the color factors 



We first consider the color operators which can appear in the quark-quark diagrams of 
sec. [IVl They read 



Oi-ijki — -^(X a )ji(X a )ik = y ■ y, 



Os-ijki — ^(A a )ij(A a )fc/ = 



- -(X a )ki(K)ji = ~ + -y 



. Ai 




aY 




2 




2 


> 


K 






2 






2 



(87a) 
(87b) 

(87c) 

(87d) 
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where the index in \ labels the two interacting particles. The first three operators are 
present in the scattering diagram of two quarks, a quark and an antiquark, and two anti- 
quarks respectively. The fourth color operator appears in the quark-antiquark annihilation 
diagram. It is logically a projector on the octet state, since this annihilation can not oc- 
cur if the quark and the antiquark are in a singlet. The color factors corresponding to the 
operators flH7|) are summed up in table [III 

If we denote as the generators of SU(3) in the representation corresponding to the 
particle i, we can check that the relation 

T* ■f?= 1 - [(f£) 2 - fr - T? 2 ] (88) 

is verified as expected for the scattering diagrams. T-jf are the generators for the repre- 
sentation of the bound state of 1 and 2. Let us note that, in eq. ( IHTl) . T 3 = A/2 and 
f 3 = -A*/2. 

The interactions between two gluons are treated in sec. [V] It appears that, among the 



TABLE I: Normalised CGC for two particles in the respective representations X and Y of SU(3), 
coupled so that they are in the representation Z. The indices i,j, k take values from 1 to 3, and 
o, b, c from 1 to 8. The color indices of the two particles are the upper indices. The lower one is 
the color of their coupling. 

[X,Y] Z Example CGC 
[3,3] 3 qq in qqq ^pf\ 
[3,3] 3 qq in qqq 

qq in qq ^ 
[3,3] 8 qqmqqg^(\ a r 
[3,8] 3 qgmqqg^X' 1 )^ 
[3,8] 3 qg in qqg ^(X a )/ 
M 1 gg in gg ±6 ab 
[8,8] 8s gg in ggg ±f ab c 
[8,8] 8 ° gg in ggg^fld ab c 
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TABLE II: Values of the color factors corresponding to the operators (|87p for the different couplings 
of tabled! The Cj are computed thanks to formula (|82[> in the relevant cases. 





Y] 


Ci 


c 2 


c 3 


C 4 


[3, 


3} 6 


-2/3 


X 


X 


X 


[3. 


3] 3 


X 


X 


-2/3 


X 


[3, 


3] 1 


X 


-4/3 


X 





[3, 


3] 8 


X 


1/6 


X 


1/2 



four diagrams, there are only two independent color operators 

O^-abcd = —f^afedc, (89a) 
O^abcd = f e ca febd- (89b) 

The first one appears in the scattering of two gluons, while the annihilation diagram is 
proportional to second one. The corresponding color factors are given in table IIII1 Our 
results agree with those of ref. |7j. We can observe that the contribution of the annihilation 
diagram is zero when the two gluons are in a singlet, as expected, but also when they are 
in a symmetric octet. This can be encountered in a glueball formed of three gluons, and a 
posteriori justified by observing that the effective potential for the annihilation diagram is 
proportional to 5 3 (r). Thus, its contribution is nonzero only for the s-wave, which has an 
symmetric space wave function and an antisymmetric spin wave function. In order for the 
total wave function to be symmetric (we are dealing with two identical bosons), the color 
function needs thus to be antisymmetric. That is why the annihilation can not occur in the 
symmetric octet. Let us also note that eq. (1881) can be used in the scattering diagram of 
two gluons, but in this case, (T 8a ) bc = —if a bc have to be used, as the generators of the 8 
representation. 

The last case concerns the quark-gluon interactions. Five different color operators appear 
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TABLE III: Same as table HH but for the color operators (|89p . 





c 5 


/7 

c 6 


[8.8] 1 


-3 





[8,8] 8s 


-3/2 





[8,8] 8a 


-3/2 


-3 



TABLE IV: Same as table HH but for the color operators (p?0|h 



[X,Y] Z 


C-7 Cg 


c 9 






[3,8] 3 


-3/2 4/3 


X 


-1/6 


X 


[3, 8] 5 


-3/2 x 


4/3 


X 


-1/6 



in the diagrams 



@7;bcjk 

Og-bcjk 
@lO;bcjk 
On-bcjk 



[-^)kj(-ifacb) 



(\ c X b ) kj , 

(\ b \ C )jk, 



1 

4 
1 

4 

^(A fe A c ) fci , 
-(X c X b ) jk , 



(90a) 
(90b) 
(90c) 
(90d) 
(90e) 



and the corresponding color factors are listed in table [TV] Again, C 7 can be computed with 
eq. ( BSD , m agreement with ref. |l2l |. 



X. CONCLUSION AND OUTLOOK 



We have computed in this paper the effective potentials arising from second order Feyn- 
man diagrams involving quarks, antiquarks, or gluons. Our expressions remain defined for 
the lightest particles, because the masses appearing in our potentials are seen as the con- 
stituent ones, not the bare ones. For instance, in the nucleon, the constituent mass of the 
u, d quarks are about 300 MeV. In glueballs, the typical value for the gluon constituent mass 
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is 500 MeV J9J . It is worth mentioning that a recent comparison between quark models and 
large N c expansion of QCD also supports the development of the relativistic corrections in 
powers of |39l |. 

We firstly re-obtained the well-known Fermi-Breit potential for two interacting quarks. 
Then, we computed the effective potential between two constituent gluons, and we derived an 
effective potential for the quark-gluon interaction and showed that in this case a quadrupolar 
interaction is present, due to the asymmetry of the system. Our results are valid at the order 
c~ 2 , which ensures us to have correct hyperfine terms, like spin-orbit, tensor force, . . . 

As some studies argue that a nonzero mass must be confered to the exchanged gluons, 
we showed how to include in a simple way such a mass starting from our expressions. 
Whether the exchanged gluons are massive or not clearly leads to very different potentials. 
The predictions of these two approaches could be compared, for example to the glueball 
spectra obtained in lattice QCD, in order to check their relevance. However, we stress that 
the goal of this work is not to decide whether the gluons are massive or not: We computed 
general effective potentials in a systematic way from QCD, and were led to expressions which 
remain valid for light particles. Furthermore, for the first time, we obtained a quark-gluon 
interaction potential. 

As an important check of the framework that we presented here, we can mention that we 
recently used it to compute light and heavy meson as well as two-gluon glueball spectra 35]. 
It can be seen in this last reference that the agreement between experimental data and 
our theoretical spectra is very satisfactory, and that the glueball spectrum we get is mostly 
located in the error bars of the lattice QCD one. Consequently, we think that the various 
potentials which are presented in this study can be successfully used as short-range terms 
within potential models describing not only usual hadrons, but also exotics like two or three- 
gluon glueballs, and hybrid qqg mesons, which currently deserve much interest. Motivated 
by the results of ref. 35], we think that our results can be seen as a basis for further 
computations of exotic hadron spectra, that we leave for future works. 

The authors thank Dr. Claude Semay for many fruitful and stimulating discussions. 
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APPENDIX A: FEYMAN RULES FOR QCD 



We recall in this section the Feynman rules we use. They are identical to those of ref. [38, 
p. 383] up to some constant normalisation factors, which were changed to obtained more 
suitable expressions for our calculations. For the in and out states wave functions, we have: 

Ingoing quark 1 S • = u l (p) 

Outgoing quark • £ = u 1 (p) 

Ingoing antiquark 1 < p • = v l (p) 

Outgoing antiquark • < p % — v l {p) 

Ingoing gluon qr \/\3\J\)» — 

Outgoing gluon *\/\jK/\? = £ a J(p), 



where p is the 4— momentum of the particle. As color indices, we used a = 1,...,8 for 
a particle in the adjoint representation, and i = 1,2,3 for a particle in the fundamental 
representation. 

The quark propagator and the gluon propagator in Lorentz gauge read 

\ + m) 



i 



iS. 



11 o 9 



VA^Ai = -iS ab D^{q). 



D^ u (q) is defined as 



DM 



(Al) 



where 9 is the gauge parameter, and r] = diag(H ) is the Minkowski metric, /i, v = 

0, ...,3 are used as space-time indices. As usual, ^ = j^q^, with 7 M the Dirac matrices, 
satisfying 

{r,7"} = 2ir- (A2) 

These matrices can be explicitly represented by 



7 



1 0^ 












[o -11 







(A3) 
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with cr the Pauli matrices. 

We also need the contributions of the different vertices. 



They are given by 





= -gfai 



■ ?) V + (<? - *0 V + {k- p)V a 



v,b,k 




where A a are the Gell-Mann matrices, and j a \, c are the structure constants of £77(3). Finally, 
for each matrix element, we have to add a factor (2tt) a 5 a (P — P') for the energy-momentum 
conservation. P and P' are the initial and final 4-momentum respectively. 



APPENDIX B: SOME FOURIER TRANSFORMS 



We present in this section useful Fourier transforms that are taken from ref. [18|, p. 282] 
We define the Fourier transform of a function g(q) as 

d 3 q 



{2kY 



where the bold quantities always denote a vector. 

Then, one can check that the following equalities hold 



(Bl) 



Air 



T 



T 

An(a • q)(b • q 



q 2 + m? 





q 4 


Air(a • 


q)(b-q) 


q 2 


+ m 2 


Aiii a 


• (q x p) 



q 2 + m 2 



S 3 (r) 



m > 0, 



r 

1 

2r 



a • b 



[a • r)(b • r) 



(a- V)(b- V) + -(a-b)V 2 

[a " (rXp)] r^^~' 



(B2) 
(B3) 

(B4) 

(B5) 

(B6) 
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with a and b vectors or matrix vectors. All these identities hold for massless propagators, 
but the limit m = has to be taken before applying the derivatives. Let us recall the relation 

(V 2 - m 2 ) — = -4vr5 3 (r), (B7) 
r 

which is useful in the computation of eq. flB5[) . 
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